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Abstract - The cancelation of the first even zonal harmonic coefficient J2 from the linear combina- 
tion /'^^^ of the nodes Q of LAGEOS and LAGEOS II used in the latest tests of the Lense-Thirring 
effect cannot be perfect, contrary to what assumed so far. It is so also because of the uncertainties 
in the spatial orientation of the terrestrial spin axis fe. As a consequence of above, the coefficient 
ci entering f''^^\ which is not a solve- for parameter being, instead, theoretically computed from 
the analytical expressions of the classical node precessions Slj^ due to J2, is, on average, uncertain 
at a 10~* level over multi-decadal time spans AT comparable to those used in the data analyses 
performed so far. A further ~ 20% systematic uncertainty, thus, occurs. The shift A/olt due 
to the gravitomagnetic frame-dragging on the station-spacecraft range p is numerically computed 
over AT = 15 d and AT = 1 yr. The need to look at such a directly observable quantity is 
highlighted, along with some critical remarks concerning the methodology used so far to measure 
the Lense-Thirring effect with the LAGEOS satellites. Suggestions for a different, more trustable 
and reliable approach are offered. 
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I Introduction. — In the Einstcinian general relativ- 
ity, which is a fully Lorcntz-invariant theory of gravitation, 
inatter-energy currents create an additional, magnetic-like 
component of the gravitational field [T] with respect to the 
static case. It is believed to play a relevant role in explain- 



J-jing relativistic jets ejected from active galactic nuclei 



-The gravitomagnetic field of a rotating body affects or- 
biting test particles, precessing gyroscopes, moving clocks 
and atoms, and propagating electromagnetic waves with 
a variety of phenomena Some of them have been put 
to the test more or less recently. For an overview of such 
a phenomenology in the solar system, see, e.g., [Fj. 

The Gravity Probe B (GP-B) experiment [S] officially 
came to an end, with the release of its final results [7] ac- 
cording to which the general relativistic gravitomagnetic 
gyroscope precession [8]|9] would have been measured with 
a claimed accuracy of 19%. Such a figure is greater than 
the previously expected 1% level because of a number of 
unwanted systematic errors whose proper treatment re- 
quired much additional efforts by the GP-B team [T^ . In- 
dependent analyses by different teams will be important 



in critically assessing the reliability of the results of [?]• 
This is beyond the scope of the present paper. It is now 
even more important than before to critically scrutinize 
the competing tests of the gravitomagnetic Lense-Thirring 
orbital precessions [TT] performed in the past years [H] 
with the LAGEOS and LAGEOS II SLR satellites in the 
gravitational field of the Earth, as originally proposed in 
[13] . Let us recall that the GP-B mission was a dedicated 
experiment in the terrestrial gravitational field costed US$ 
750 million and lasted 52 yr, while the gravitomagentic 
data analyses [13] of the LAGEOS spacecrafts, which were 
originally launched for different purposes, were much less 
expensive and comparatively less extended in time. Ac- 
cording to I. Ciufolini, its accuracy would be 10% or bet- 
ter [TSj; for recent articles establishing a comparison be- 
tween GP-B and the previous LAGEOS-based results, see 
[T5HT5] in which it is basically argued that GP-B would 
have just reached the same results of the earlier tests with 
the LAGEOS spacecraft at a much higher cost and with 
an even worst, or, at most, comparable, accuracy. 
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Can the cancelation of the effect of the 
quadrupole mass moment of the Earth in the 
LAGEOS-based tests be perfect?. - The follow- 
ing linear combination of the longitudes of the ascending 
nodcflll Q of LAGEOS and LAGEOS II pln2D] 



/(2L) = |2{L)+c^^(L II) 



(1) 



was used in the tests of the Lense-Thirring effect per- 
formed so far with such artificial bodies orbiting the Earth. 
Frame-dragging was purposely not modeled [TS] , and time 
series o^ "residuals" of the nodes [3TJ[35] of both satel- 
lites, combined according to eq. ([T]), were analyzed and 
subsequently fitted with a straight line plus other time- 
dependent signals. The coefficient c\ entering eq. ([T]) is 
not one of the several solve-for parameters estimated in 
the data reduction process. Following an approach set 
forth in a different context [52, its value is theoretically 
computed as [20] 



Cl 



(2) 



from the analytical expressions of the classical secular 
node precessions Qj^ of both the LAGEOS sateUites 
caused by the first even zonal harmonic coefficient Ji of 
the expansion in multipoles of the Newtonian part J7n of 
the terrestrial gravitational potential. This multipolar ex- 
pansion of accounts for its departure from spherical 
symmetry because of the centrifugal deformation due to 
the Earth's diurnal rotation [53]. Traditionally, eq. © 
has always been computed [50] so far from the well known 
expression 
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3nJ2R^ cos/ 
2a2(l -e2)2' 



(3) 



In eq. ([3]) a is the semi-major axis of the satellite's orbit, 
e is its eccentricity, / is its inclination to the reference 
{X, Y} plane, assumed to be coincident with the Earth's 
equator, R is the terrestrial equatorial radius, and n = 
y/ GM/a^ is the Keplerian mean motion of the satellite 
with respect to the Earth whose mass is denoted with M\ 
G is the Newtonian constant of gravitation. The orbital 
parameters of LAGEOS and LAGEOS II, referred to a 
geocentric inertial system, are shown in table [1] 

The aim of eq. ([1]), with eq. ([5]), is to cancel out, by 
construction, such precessions. Since they are nominally 7 
orders of magnitude greater than the Lense-Thirring ones 
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c2a3(l - e2)3/2' 



(4) 



^ The longitude of the ascending node Q is one of the angles de- 
termining the orientation in space of the satellite's Keplerian ellipse. 

^The term "residual" is, actually, improper for the node. In- 
deed, all the Keplerian orbital elements are not observable quanti- 
ties. They can only be computed at various epochs from the corre- 
sponding state vectors in cartesian coordinates which, in turn, are 
computed from the measured values of the direct observables. 



where S is the Earth's angular momentum and c is the 
speed of light in vacuum, they represent a major source 
of systematic bias in determining them. For the LAGEOS 
satellites eq. ([4]) yields about 30 milliarcseconds per year 
(mas yr~^ in the following), so that the combined Lense- 
Thirring signal amounts to approximately 50 mas yr~^ ac- 
cording to eq. ([IJ. In principle, such a removal of J2 from 
eq. ([U is exact, or so it has always been considered until 
now. Indeed, in all the more or less realistic evaluations 
of the systematic error due to the gcopotential existing in 
literature [51[T51[T3], the part due to J2 was always set to 
zero by definition and independently of aj^ . Thus, the fo- 
cus was on the impact of the other, uncanceled even zonal 
harmonics of higher degree Ji, ^ = 4, 6, 8, . . ., known with 
a certain level of uncertainty. Actually, the effect of J2 
on eq. ([IJ cannot be exactly zero because of a number of 
factors. 

One of them relies on the fact that, for a given set of 
valued of the satellites' orbital parameters from which ci 
is computed by means of eq. ([2]) and eq. ([3]), the ac- 
tual accuracy with which ci can be known is necessarily 
limited by the uncertainties with which the satellites' Ke- 
plerian orbital elements of interest can be determined in 
the data reduction procedure. It was recently shown |25j 
that CTa — 2 cm and u/ ~ 0.5 mas yield Aci ~ 10~*, cor- 
responding to a further systematic uncertainty of about 
20% in the Lense-Thirring signature. If, instead, one op- 
timistically assumes Ca ~ 2 cm and [12] c/ ~ 10 — 30 //as, 
then Aci ~ 8 x 10^^, which implies an additional 14% 
bias. On the contrary, ci was always released so far with 
a very limited number of significant digits; for example, in 
[T5] we have ci = 0.545. As pointed out in [25^, it would 
be incorrect to argue that the impact of Aci would be 
negligible since it should be multiplied by the uncertainty 
in J2. Indeed, the standard error propagation theory tells 
us that, in addition to the mixed, cross-correlated terms 
containing the products of the uncertainties, there are also 
the linear terms proportional to the uncertainties in each 
parameter. Moreover, in the LAGEOS tests both ci and 
J2 are not estimated solve-for parameters. For the sake of 
definitcness, we will denote the values of ci obtained from 
eq. © with cf^; tabled] yields c^i^ = 0.540976405. 

Another issue, not yet considered in literature, is that it 
is incorrect to assume a perfect alignment of the Earth's 
spin axis, whose unit vector is denoted by k, and the refer- 
ence Z axis of the geocentric inertial system actually used. 
Indeed, on the one hand, the latter refers to a given ref- 
erence epoch, typically J2000.0, while the time spans AT 
over which the data of LAGEOS and LAGEOS II were an- 
alyzed necessarily cover 19 yr or less: during such a tempo- 
ral interval k did not remain fixed in the inertial space due 
to a variety of physical processes [5B]. Such changes, even 
if taken into account and modeled, are, of course, known 



•^They were never explicitly specified in the analyses performed so 
far, by assuming for them some standard figures 1121 . approximately 
representative of the orbital configurations of the LAGEOS satellites. 
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Table 1; Keplerian orbital parameters of LAGEOS and LAGEOS II computed from state vectors, in cartesian coordinates, 
corresponding to a given epoch kindly provided by L. Combrinck to the author. The inclination / and the node Q refer to a 
geocentric inertial system whose reference {X, Y} plane is assumed to be coincident with the Earth's equator. The semi-major 
axes and the angles are given with a cm-level and mas-level accuracy, respectively (1 cm= 10^'' km, 1 mas= 2.7 x 10^^ deg). 

Spacecraft a (km) e / (dcg) Q (deg) 

LAGEOS 12274.75303 0.0039962 109.8617388 -1.4477848 
LAGEOS II 12159.19724 0.0141892 52.6013013 -94.7543331 



only with a limited accuracy [27] . On the other hand, it is 
well known that another source of uncertainty in the loca- 
tion of fc is given by the polar motion |26| with respect to 
the Earth's crust itself, known with an accuracy of about 
10 — 20 mas [5S1I2H] over a time interval of just 1 yr. See 



Uncertainty in cj (fix = Sy = 10 mas) 



also http://www.iers.org/nn_10398/IERS/EN/Science/ 



EarthRotation/PolarMotion.html?_nnn=true. Thus, it is 
important to quantitatively assess the further systematic 
error Aci induced by the use of c^^ with respect to val- 
ues, denoted as c'i'°\ computed by taking into account the 
real spatial orientation of k. To this aim, a first step con- 
sists of computing the long-term node variations Qj^ for a 
generic orientation of k. The acceleration experienced by 
a test body orbiting an oblate central mass rotating about 
a generic direction k is pS] 



lJ2 



3GMJ2R^ 
2^ 



1 - 5 f • fc 



2 r • fc fe 



(5) 

Since its magnitude is quite smaller than the main Newto- 
nian monopole, its effect on the particle's orbital motion 
can be straightforwardly worked out with standard pertur- 
bative techniques. The Gauss equation for the variation 
of the node (SS] allows to obtain the rate of change of H 
averaged over one orbital revolution. It turns out to be 

- . -F(/,/?;fc), (6) 



with 



4a2 (1 - e^y 



2kz cos 21 CSC / ( kx sin Q — ky cos 17 
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k\ 



k\ \ COS 2/2- 



2kxkY sin2J7 



/7) 

It is an exact result in e and / in the sense that no a-priori 
simplifying assumptions on their values were assumed; in 
general, it can also be useful in other contexts involving 
different central bodies and test particles [S] . It can be no- 
ticed that, according to eq. © and eq. ([7]), the long-term 
rate of change of Q consists of the sum of a genuine sec- 
ular precession and of a harmonic, time-dependent signal 
involving fl and 2Q. Moreover, eq. ([7]) reduces to 

(8) 
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Fig. 1: Temporal evolution of the difference Aci between the 
value of cj"' computed by assuming the Earth's spin axis ex- 
actly coincident with the reference Z axis of an inertial equa- 
torial reference system, and the value of c-^^'^ computed by 
assuming an uncertainty of the order of 10 mas in the ori- 
entation fc of the Earth's spin axis in the same reference system. 
The time span is AT = 19 yr, while the time step is At = 7 
d. The initial conditions chosen for LAGEOS and LAGEOS 
II are those listed in table [T] Aci is characterized by an av- 
erage (Aci) = 4.1 X 10^**, and a peak-to-peak amplitude of 



for kx ~ ky = 0,kz = ±1, yielding the well-known sec- 
ular precession of cq. ([3]) . We will denote the value of ci 



computed from eq. ©-eq. ([7]) by 



In Figure [T] we 



plot the uncertainty in ci raising from having used just 
cf'^ over a temporal interval AT = 19 yr representative 
of the time spans actually used in real data analyses, and 
for a 10 mas uncertainty in the position of fc. It can be 
noticed that its impact is non-negligible since it is of the 
order of 4 — 8 x 10^^, implying a further ~ 20% systematic 
uncertainty in the gravitomagnetic signature. 

What was really measured in the LAGEOS- 
based tests?. — In SLR studies, the directly observ- 
able quantity is the range p between a spacecraft equipped 
with retroreflectors and a ground-based statiorQ . It is 
straightforwardly computed by multiplying c by the time 
interval elapsed between the emission of the laser pulse 



T = -2 cos/ 



*It is just the case to recall that ranges refer to an Earth-fixed 
rotating reference system. In order to obtain the values of table ^ 
one has to take into account the polar motion, the Earth rotation, 
the precession and the nutation. See I26| for details. 
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sent to the orbiting target body and its subsequent recep- 
tion after it was bounced back by the retroreflectors on- 
board the sateUite. The precision of such measurements 
is nowadays at the mm level [26|. Post- fit range residuals 
for good targets like LAGEOS and LAGEOS II, obtained 
after the adjustment of a number of solvcd-for parameters 
pertaining to the satellites' physical properties and orbital 
dynamics, and the measurement process itself, are as large 
as 1 cm or less in a Root-Mcan-Square (RMS) sense [26] . 
They globally reflect the impact of all the unmodclcd and 
mismodeled sources of errors like, e.g., some unknown or 
poorly modeled forces acting on the satellites. The post- 
fit range residuals are also a measure of the effectiveness 
of the orbit determination process in which the estimated 
values of some parameters may partly or totally absorb the 
effects of other parameters not included in the list of those 
to be adjusted, or of totally unmodclcd forces themselves. 
In general, if one is interested in a certain dynamical fea- 
ture, then it must be explicitly modeled in such a way 
that one or more dedicated solvc-for parameters are esti- 
mated. Subsequently, the resulting covariance matrix can 
be examined to identify the correlations between various 
parameters. Clearly, the magnitude of post-fit range resid- 
uals can only be greater than, or as large as the range mea- 
surement precision. Perfect models and/or total removal 
of all effects that have not been modeled would provide 
residuals as large as the measurement precision. 

Extending such considerations to the frame-dragging 
tests made so far with the LAGEOS satellites, it must 
be remarked that, actually, the Lense-Thirring force was 
never modeled, so that it should be considered in the same 
way as a source of systematic error impacting, in case, the 
post-fit range residuals to a certain level. No dedicated 
solve-for parameters were ever estimated; thus, the grav- 
itomagnetic signature might have been partly or totally 
absorbed in the estimation of the several other parame- 
ters in the data reduction process, and partially or to- 
tally removed from the range signature. If frame-dragging 
fully impacted the ranges as predicted by general relativ- 
ity, there should be time series of post-fit range residuals 
with the characteristic signature of the gravitomagnctic 
force itself. See Figure [2] and Figure [3] displaying the nu- 
merically produced nominal Lense-Thirring effect on the 
station-satellite range for LAGEOS and LAGEOS II over 
a time interval of AT = 1 yr. On the other hand, the same 
set of data should be analyzed by explicitly modeling the 
Lense-Thirring effect in order to check if statistically sig- 
nificant differences with respect to the previous case would 
occur. This would be a crucial test of the ability to ac- 
tually measure terrestrial gravitomagentism by means of 
the LAGEOS and LAGEOS II SLR data. In fact, after 
more than 15 years since the first tests, such "gravito- 
magnctic" post-fit range residuals were never shown so 
far. It should be noticed that there is a contradiction be- 
tween claiming sub-cm post-fit range-residuals, obtained 
without modeling frame-dragging, and figs. [2] and [3] dis- 
playing signatures with RMS variances as large as 18.0 



Station-LAGEOS Lense-Thirring range sliift 
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Fig. 2: Numerically integrated Lense-Thirring station-satellite 
range perturbation ApLT for LAGEOS over AT — 1 yr. Its 
variance is 18.0 cm. We choose the ITRF2000 coordinates of 
the GRAZ station, from [30]: cut-off elevation angle of 20 deg. 



cm and 46.1 cm, respectively. Indeed, one should assume 
either that the gravitomagnctic signal, not modeled, was 
almost entirely removed or that it was almost canceled by 
the superposition of other unmodclcd/mismodcled com- 
peting dynamical effects. After all, such a removal would 
not be implausible since, as shown by fig. H] and fig. [5j 
the nominal size of the Lense-Thirring range perturbation 
is just at the level of cm on a timescale of At = 15 d. It 
is not clear, however, why all the other effects not mod- 
eled at all, or poorly modeled, should be exactly removed, 
or should cancel each other leaving just the completely 
unmodclcd Lense-Thirring signal, which is precisely what 
one expects to find in the data. It is much more plau- 
sible that it is somewhat absorbed in some of the esti- 
mated parameters and removed from the residual signal 
to a certain extent. Somebody may argue that the re- 
moval of the Lense-Thirring signature can occur only if 
certain once-per-revolution empirical cross-track acceler- 
ations were estimated. First of all, it should be explic- 
itly proven that they were actually not estimated in the 
dedicated LAGEOS data reductions. More importantly, 
it is impossible to a-priori decide in which of the esti- 
mated parameters the cancelation would actually occur. 
Suffice it to say that in much more "clean" scenarios like 
planetary astronomy, not plagued by the host of distur- 
bances and non-gravitational effects of satellite geodesy, 
it is common practice to explicitly model the effects one 
is interested in and solve for one or more dedicated pa- 
rameters just to avoid the risk that they may be partially 
or totally absorbed in the estimation of the initial state 
vectors. Interestingly, this has been done recently [5TH5^ 
even for hypothetical forces that, as the Pioneer Anomaly, 
if they really existed in Nature would have caused signa- 
tures much greater than the accuracy of the observations 
themselves. 

We remark that the LAGEOS-based tests are likely 
plagued by another source of intrinsic a-priori imprint- 
ing of general relativity itself in addition to those already 
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Station-LAGEOS II Lense-Thirring range shift 
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Fig. 3: Numerically integrated Lense-Thirring station-satellite 
range perturbation ApLT for LAGEOS II over AT = 1 yr. Its 
variance is 46.1 cm. We choose the ITRF2000 coordinates of 
the GRAZ station, from 'SCT: cut-off elevation angle of 20 deg. 



Station-LAGEOS Lense-Thirring range shift 
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Fig. 4: Numerically integrated Lense-Thirring station-satellite 
range perturbation ApLx for LAGEOS over AT — 15 d. Its 
variance is 0.8 cm. We choose the ITRF2000 coordinates of 
the GRAZ station, from [30]: cut-off elevation angle of 20 deg. 



Station-LAGEOS II Lense-Thirring range shift 
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Fig. 5: Numerically integrated Lense-Thirring station-satellite 
range perturbation ApLT for LAGEOS II over AT — 15 d. Its 
variance is 1.9 cm. We choose the ITRF2000 coordinates of 
the GRAZ station, from ,30, : cut-off elevation angle of 20 deg. 



pointed out [5]- Indeed, they always made use of a refer- 
ence system whose materialization heavily relies upon SLR 
data, among which those from LAGEOS and LAGEOS II 
themselves play a fundamental role. 

The considerations exposed here are, in principle, valid 
also for other performed or proposed tests of general rela- 
tivity with the LAGEOS satellites [3S], and also for those 
which should be implemented in the near future with the 
existing LAGEOS and LAGEOS II, and with the new 
LARES satellite [36], to be launched in late 2011 with 
a VEGA rocket. 



Conclusions. — In conclusion, we can entertain rea- 
sonable doubts as to what it was actually seen in the tests 
with LAGEOS and LAGEOS II made so far, and what 
has been passed of as frame-dragging in them. Only the 
use of a completely different approach, more related to 
quantities that are actually measured, could afford to talk 
about of clear and unambiguous tests of this subtle ef- 
fect. Frame-dragging should be explicitly modeled and 
solved-for in the LAGEOS and LAGEOS II data reduc- 
tion process; post-fit range residuals produced with and 
without a model for the Lense-Thirring effect should be 
displayed and analyzed; a different materialization of the 
reference system used so far, mostly based on the observa- 
tions of LAGEOS and LAGEOS II themselves, should be 
adopted; it would be preferable that GR is explicitly mod- 
eled and solved-for in future dedicated global gravity field 
solutions combining data from several satellites. Other- 
wise, they should make clear why they do not implement 
the strategy advocated here which, after all, is standard 
practice in all branches of geodetic and astronomical stud- 
ies. 

Moreover, even accepting the strategy followed so far, 
the unavoidable uncertainties in our knowledge of the 
Earth's rotation axis affect the necessarily imperfect cal- 
culation of the theoretical coefficient ci entering the linear 
combination of the nodes of LAGEOS and LAGEOS II. 
It does not allow to obtain an exact cancelation of the 
aliasing bias due to the first even zonal harmonic J2 of 
the geopotential which, instead, would still be present at 
a ~ 20% level of the Lense-Thirring signal. Let us recall 
that a further 10 — 20% alias comes from the uncertainty 
in ci due to the errors in the satellites' orbital parameters 
a and /. 

Thus, more work is still needed to really consider the 
LAGEOS-based attempt as a robust complement of the 
GP-B mission from the point of view of reliability, trusta- 
bility and methodology. Although the LAGEOS-based 
tests had measured something that really relates to the 
Lense-Thirring effect, their overall uncertainty will prob- 
ably make them less accurate than the GP-B experiment. 
Anyway, independent analyses of the data of the Stanford 
team by different groups are certainly required. 
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* * * 

I thank L. Combrinck for having provided me with the 
state vectors of both the LAGEOS spacecrafts, and for 
interesting discussions. 

REFERENCES 

[I] Thorne K. S., MacDonald D. A. and Price R. H. (Edi- 
tors) , Black Holes: The Membrane Paradigm (Yale Univer- 
sity Press, Yale) 1986. 

[2] Thorne K. S., Near Zero: New Frontiers of Physics, 
edited by J. D. Fairbank, B. S. Deaver, C. W. F. 
Everitt and P. F. Michelson (W. H. Freeman and Com- 
pany, New York) 1988, p. 573. 

[3] Stella L. and Possenti A., Space Sci. Rev., 148 (2009) 
105. 

[4] SCHAFER G., Space Sci. Rev., 148 (2009) 37. 

[5] lORIO L., LiCHTENEGGER H. I. M., RUGGIERO M. L. and 

CORDA C, Astrophys. Space Sci., 331 (2011) 351. 
[6] Everitt C. W. F., Buchman S., DeBra D. B., et 

al, Gyros, Clocks, Interferometers...: Testing Relativis- 

tic Gravity m Space, Vol. LNP 562, edited by C. 

Lammerzahl, C. W. F. Everitt and F. W. Hehl 

(Springer, Berlin) 2001, p. 52. 
[7] Everitt C. W. F., DeBra D. B., Parkinson B. 

W., TURNEAURE J. P., CONKLIN J. W., HeIFETZ M. 

I., Keiser G. M., Silbergleit A. S., Holmes T., 
KoLODZiEJCZAK J., Al-Meshari M., Mester J. C., 

MUHLFELDER B., SOLOMONIK V., StAHL K., WORDEN 

p., Bencze W., Buchman S., Clarke B., Al-Jadaan 
A., Al-Jibreen H., Li J., Lira J. A., Lockhart J. M., 
Al-Suwaidan B., Taber M., and Wang S., Phys. Rev. 
Lett, 106 (2011) 221101. 

[8] Pugh G. E., WSEG Research Memorandum, Vol. 11 (The 
Pentagon, Washington 25 D.C.) 1959. 

[9] Schiff L. I., Phys. Rev. Lett., 4 (1960) 215. 

[10] Muhlfelder B., Adams M., Clarke B., Keiser G.M., 

KOLODZIEJCZAK J., Ll J., LOCKHART J. M., WORDEN P., 

Space Sci. Rev., 148 (2009) 429. 

[II] Lense J. and Thirring H., Phys. Z, 19 (1918) 156. 
[12] CiUFOLiNi I., Paolozzi A., Pavlis E. C, Ries J. C, 

KONIG R., Matzner R. a., Sindoni G. and Neumayer 
H., Space Set. Rev., 148 (2009) 71. 
[13] CuGUSi L. and Proverbio E., Astron. Astrophys., 69 
(1978) 321. 

[14] CiUFOLiNi I., Pavlis E. C, Ries J. C, Konig, R., 
Sindoni G., Paolozzi A. and Neumayer H., General 
Relativity and John Archibald Wheeler, Vol. ASSL 367, 
edited by I. Ciufolini and R. A. Matzner (Springer, 
Dordrecht) 2010, p. 371. 

[15] Ciufolini I. and Pavlis E. C, Nature, 431 (2004) 958. 

[16] Reich E. S., Nature, 473 (2011) 131. 

[17] Cho a.. Science NOW, (2011) 

http:/ /news.sciencemag.org/sciencenow/2011/05/at- 
long-last-gr avity- probe-b.html. 

[18] Reich E. S., Nature News Blog, (2011) 
http:/ /blogs. nature.com/news/2011/05/post_76. html. 

[19] Ries J. C, Eanes R. J., Tapley B. D. and Peter- 
son G. E., Proceedings of the 13th International Laser 
Ranging Workshop, edited by R. Noomen, S. Klosko, 



C. Noll and M. Pearlman, Vol. NASA CP (2003- 
212248) (NASA, Goddard) 2003. 

[20] lORiO L., THE TENTH MARCEL GROSSMANN 
MEETING. On Recent Developments m Theoretical and 
Experimental General Relativity, Gravitation and Rela- 
tivistic Field Theories. Proceedings of the MGIO Meeting 
held at Brazilian Center for Research in Physics (CBPF), 
Rio de Janeiro, Brazil, 20-26 July 2003, edited by M. Nov- 
ELLO, S. P. Bergliaffa and R. J. Ruffini (World Scien- 
tific, Singapore) 2006, p. 1011. 

[21] LUCCHESI D. M. and Balmino G., Planet. Space Sci., 54 

(2006) 581. 

[22] LucCHESi D. M., Adv. Space Res., 39 (2007) 1559. 
[23] Ciufolini I., Nuovo Cimento A, 12 (1996) 1709. 
[24] Capderou M., Satellites. Orbits and missions (Springer, 
Paris) 2005. 

[25] lORlO L., Gen. Relativ. Gravit., 43 (2011) 1697. 

[26] Montenbruck O. and Gill E., Satellite Orbits. Models, 
Methods, and Applications (Springer, Berlin) 2000. 

[27] Seidelmann P.K., Archinal B.A., A'Hearn M.F., 
Conrad A., Consolmagno G.J., Hestroffer D., 
Hilton J.L., Krasinsky G.A., Neumann G., Oberst 
J., Stooke p., Tedesco E.F., Tholen D.J., Thomas 
P.C., and Williams LP., Celest. Mech. Dyn. Astron., 98 

(2007) 155. 

[28] McCarthy D. D. and Petit G., IERS Conventions 
(2003) (Verlag des Bundesamtes fiir Kartographie und 
Geodasie, Frankfurt am Main) 2004. 

[29] Bertotti B., Farinella P. and Vokrouhlicky D., 
Physics of the Solar System (Kluwer, Dordrecht) 2003. 

[30] Zhao G., Zhao Y., Sun M., Yu H., 16th In- 
ternational Workshop on Laser Ranging, (2008) 
http: / /cddis. gsfc.nasa.gov /lwl6 / docs /presentations / net_9_Zhao.pd 

[31] Standish E. M., AIP Conference Proceedings: Recent 
Developments in Gravitation and Cosmology-3rd Mexican 
Meeting on Mathematical and Experimental Physics, Vol. 
977, edited by A. Macias, C. Lammerzahl, and A. Ca- 
MACHO (American Institute of Physics) 2008, p. 254. 

[32] Standish E. M., Relativity in Fundamental Astronomy: 
Dynamics, Reference Frames, and Data Analysis, Proceed- 
ings lAU Symposium No. 261, edited by S. A. Klioner, 
P. K. Seidelmann, and M. H. Soffel (Cambridge Uni- 
versity Press, Cambridge) 2010, p. 179. 

[33] Fienga a., Laskar J., Kuchynka P., Leponcin- 
Lafitte C, Manche H., and Gastineau M., Relativity 
in Fundamental Astronomy: Dynamics, Reference Frames, 
and Data Analysis, Proceedings lAU Symposium No. 261, 
edited by S. A. Klioner, P. K. Seidelmann, and M. H. 
Soffel (Cambridge University Press, Cambridge) 2010, 
p. 159. 

[34] FOLKNER W. M., Relativity in Fundamental Astronomy: 
Dynamics, Reference Frames, and Data Analysis, Proceed- 
ings lAU Symposium No. 261, edited by S. A. Klioner, 
P. K. Seidelmann, and M. H. Soffel (Cambridge Uni- 
versity Press, Cambridge) 2010, p. 155. 

[35] LucCHESi D. M. and Peron R., Phys. Rev. Lett., 105 
(2011) 231103. 

[36] Ciufolini L, Paolozzi A., Pavlis E. C, Ries J. C, 
KoNiG, R., Matzner R. A. and Sindoni G., General 
Relativity and John Archibald Wheeler, Vol. ASSL 367, 
edited by I. Ciufolini and R. A. Matzner (Springer, 
Dordrecht) 2010, p. 467. 



p-6 



